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Abstract 

Discrete multiparametric ID quantum well with VT —symmetric long-range boundary condi¬ 
tions is proposed and studied. As a nonlocal descendant of the square well families endowed 
with Dirac (i.e., Hermitian) and with complex Robin (i.e., non-Hermitian but still local) bound¬ 
ary conditions, the model is shown characterized by the survival of solvability in combination 
with an enhanced spectral-design flexibility. The solvability incorporates also the feasibility 
of closed-form constructions of the physical Hilbert-space inner products rendering the time- 
evolution unitary. 
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1 Introduction 


An active interest of researchers in quantum lattices alias quantum chain models alias tridiagonal- 
matrix Hamiltonians 
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dates back to the pioneering applications of these Hamiltonians in organic chemistry in the 
thirties [Tj. Still, their study belongs to the mainstream activities, say, in the tight-binding 
descriptions of electronic structures in solids [2], etc. The models of this type also helped to 
clarify some non-variational features of Green’s functions in particle physics [3] and they are 
currently serving as an exemplification of several gain-and-loss-related phenomena in optics p[j. 
Last but not least, the use of models (JTJ) threw new light on some questions of the emergence 
and confluence of the Kato’s exceptional points (KEP, [5j) in perturbation theory |6j[7]. 

In our recent paper [8] we decided to pay a detailed attention to one of the simplest quantum 
models of this type, characterized by the next-to-trivial (N +1) by (N +1) matrix Hamiltonian 
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with a single complex non-real parameter z e C \ M. We revealed that in spite of its manifest 
non-Hermiticity the model may be perceived as belonging to the conventional quantum theory 
in which it generates a unitary evolution of the system in question. 

The results of paper [8] offered an explicit example and illustration of certain recent, 
not entirely conventional implementations of the abstract principles of quantum theory. In 
essence, these implementations are based on an innovative parity-times-time-reversal-symmetric 
(VT -symmetric, PTS, [9]) presentation alias pseudo-Hermitian representation (PHR, [TO] ) 
alias three-Hilbert-space (THS, |llj, T2]) form of the formalism of quantum theory (see Ap¬ 
pendix A for a compact summary of these ideas). 
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In paper [8j we just confirmed that in practice, the costs of the PTS/PHR/THS enhancement 
of the flexibility of models may be reasonable and acceptable. Our main result was that 
we managed to fulfill the technically most difficult task of the PTS/PHR/THS theory and 
reconstructed the nontrivial “standard” representation of the physical Hilbert space in 
which the manifestly non-Hermitian Hamiltonian (J2J) becomes Hermitized so that, in other 
words, the evolution in time acquires the standard unitary interpretation. 

A formal simplicity of the discrete model (|2]) may be perceived as inherited from its 
differential-operator predecessor of Ref. [IB]. In parallel, it should be emphasized that these 
two models also share the physical motivation and an immediate phenomenological appeal 
which was verbalized in Ref. [IT] and which appeared to lie in the existence of a close relation¬ 
ship between the bound-state and scattering experimental data (cf. also Refs. [Id, 16] in this 
context). 

All these observations provoked, naturally, a search for the generalizations which would 
eventually go beyond the tridiagonal matrix structure ([Tj) of the model. A few new results 
obtained in this direction will be presented in what follows. First of all, our inspiration by PTS 
matrix ([2]) will lead, in a way described in section [2] to its partitioning and to its subsequent 
multiparametric PTS generalization. In section [3] we shall illustrate some descriptive merits of 
spectra provided by this generalization. The survival of solvability of the model will finally be 
demonstrated, constructively, in section [4j 

Due to the non-Hermiticity of the Hamiltonian, we will have to address also the above- 
mentioned problem of construction of the correct physical representation Hilbert space. This 
will be done in several sections. Firstly, the general recipe will be presented in section [5] Two 
alternative methods of construction of the correct inner-product metric 0 will be described. 
We shall point out that both parts of the construction of the bound state solutions (viz., 
the construction of the wave functions and of the metric) appear closely interconnected in a 
way which may be perceived as an explanation why, in spite of its enhanced flexibility, our 
generalized model remains tractable non-numerically. 

An alternative approach capable of producing some of the metrics in a friendlier sparse- 
matrix form will be then conjectured and studied in sections [6] and [7] In the former section we 
shall put more emphasis upon the phenomenological aspects of the metrics reflected by specific 
choices of optional parameters. In the latter section we point out that the advantages of the 
sparse metrics are accompanied by some obstacles. They will be exemplified by the emergence 
of certain critical parametric hypersurfaces of the loss of positivity of the metric, i.e., in the 
language of physics, of a phase-transition loss of the usual probabilistic interpretation of the 
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model in question. 

In section [8] we then add a brief remark on what happens, beyond the realm of quantum 
theory, when the invertibility survives but the positivity of the metric 0 is lost. Purely formally, 
matrix 0 then acquires a new meaning of pseudometric V. Although it still renders the given 
Hamiltonian self-adjoint in an ad hoc Krein or Pontriagin space ra, its physical meaning can 
only survive, say, in some less ambitious, non-unitary effective-theory implementations [T8I119J. 

The last section [9] will offer a brief summary of our results. In order to keep our present paper 
sufficiently self-contained we decided to complement the main text by two Appendices. In the 
first one we recall Ref. [T2jJ and summarize the key features of the PTS/PHR/THS quantum 
mechanics as a formalism which is based on the simultaneous use of a triplet of alternative 
representation Hilbert spaces {T-L (P \ 'Hd h \ "H^}. In Appendix B we shall turn attention to 
the well known (though rarely emphasized) specific fact that in an experimental setting, any 
generic PTS/PHR/THS model controlled by a non-Hermitian Hamiltonian is in fact nonlocal. 
This observation may be read as an idea which, initially, motivated our present considerations. 


2 The model 


2.1 Local boundary conditions 

In virtually any textbook on conventional quantum mechanics one finds the formalism illus¬ 
trated by a one-dimensional square well. In such a model the wave functions obey differential 
Schrodinger equation 

d 2 

- = E n il) n {x) (3) 

and the Dirichlet boundary conditions are imposed at a pair of boundary points, 


ij) n (±L) = 0 , L e ffi \ {0} . 


(4) 


hi pSj, similarly, the PTS/PHR/THS formalism was illustrated using a replacement of Eq. (J4]) 
by complex Robin boundary conditions 

H±L) = * -£-i/j(±L) , a,/3eR. (5) 

a =F i p dx 

In this context our subsequent paper [8] merely replaced differential Eq. (J3]) by its difference- 
equation alternative 
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We emphasized there that one of benefits of the discretization (l3li —> (J6j) lies in a survival of 
simplicity of boundary conditions in which two Eqs. (J5]) were merely replaced by their two 
equally elementary discrete analogues 

= .#& = ** AT m 

z = z(i, 0 = (i-5-i0- 1 eC, «,CeR. 

The role of the original continuous real coordinate x G K. was smoothly transferred to its 
discrete real counterpart k e Z. As a consequence the Hamiltonian itself degenerated to the 
above-mentioned finite-dimensional matrix Hamiltonian (]2j) . 

2.2 Long-range boundary conditions 

Hamiltonian (J2]) has the conventional form of superposition H = T + V . In contrast to its 
standard kinetic-energy part T, the form of the local (i.e., diagonal) potential-energy matrix 
V is much less standard. Not only that it represents just a point-supported interaction (i.e., 
it vanishes everywhere except the end-points of the discrete-coordinate lattice) but it is also 
manifestly non-Hermitian, V ^ V'. 

In papers [I8j [2D] it has been shown that in the ultimate experimental setup (i.e., in the 
primary physical Hilbert space both the locality and the point-interaction nature get 

inadvertently lost for a generic non-Hermitian interaction V, local or nonlocal, (in detail the 
reasons are recalled in Appendix B below). In other words, the decisive reason for the preferred 
choice of matrices V shouldn’t be seen in their diagonality (i.e., in their locality without physical 
meaning - cf. comments in [18]) but merely in their sparsity. Guaranteeing that with the growth 
of the matrix dimension (N + 1), the number M(V) of non-vanishing matrix elements becomes 
negligible in comparison with the total number (N + l) 2 of all matrix elements. 

Naturally, the latter condition is more than perfectly satisfied by the M(V ) = 2 Hamiltonian 
(J2D of paper PJ. We believe, nevertheless, that the restriction M(V) = 2 may be weakened at 
a reasonable cost. In our present sequel of the latter study we intend to generalize boundary 
conditions to a nonlocal constraint, therefore. As long as we still intend to keep contact with 
their local predecessors, we shall postulate that if I 7^ ^ 0 then either % or j is equal to 0 or N. 

The latter postulate admits all complex interaction matrices of the following partitioned 
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form 



bo 

* 

«i 

* 

a 2 . 

* 

■ a N -1 

do 


h 

0 

0 . 

0 

di 

u (JV+1) (d,fe,c,d) = 

b 2 

0 

0 . 

0 

d 2 


&JV-1 

0 

0 . 

0 

dN- 1 


b N 


<3 ■ 

■ C iV-l 

d N 


Obviously, this yields a class of multiparametric sparse-matrix models with M(V) = O(N). In 
parallel, in a less formal and more phenomenologically oriented perspective a truly remarkable 
feature of such forces may be seen in their surface-related support leaving the internal, bulk 
part of our discrete lattice uninfluenced. In this sense, models (jHJ) may be perceived as a certain 
missing link between the strictly local and some more general and realistic, spatially smeared 
thick-boundary-layer interactions. 


2.3 VT— symmetry constraint 


In conventional quantum mechanics one would further restrict the set of parameters in Eq. ([8]) 
by the Hermiticity constraint V — V"k This means that with a = bo, Cg = b]y, a* N = do 
and c* N = d N one would impose the two additional (N + 1)—plets of requirements a = b and 
c = d. In our present paper we shall proceed in different, PTS/PHR/THS spirit. In the 
light of multiple related technical challenges we shall, first of all, try to support the reality 
of the spectrum of H ^ H 1 by an additional, auxiliary PTS condition HVT = VTH alias 
parity-pseudo-Hermiticity constraint II'V = VH. 

Using the antidiagonal parity-representing matrix 


-p — p>(N+l) 
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01 0 ... 0 
1 0 ... 0 0 
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this leads to the most general PTS Hamiltonian 





( 10 ) 


In the notation of Eq. (J8]) this matrix contains the real elements b^ = cl and do = b and 
the mutually conjugate complex pair of 6 0 = — z and djy = —z*. The new symbols a and f3 
now denote the (N — 1)—dimensional complex vectors such that the M(V) = 2 non-Hermitian 
model 02D °f paper [8] re-emerges here as the simplest special case with a = b = 0 and a — 0 
and = 0. 


3 Spectral-design applications 

In contrast to the matrix-boundary models (JS]) which are chosen Hermitian and which would 
generate the robustly real spectra, the violation of the Hermiticity may render the parametric 
dependence of the spectrum much more versatile. Typically, for our present VT —symmetric 
model (HU the reality of the levels may be lost at the KEP values of parameters. Thus, in 
principle, the model may be expected to offer a wealth of unusual spectral patterns. Also the 
multiparametric nature of our matrix-boundary conditions would enhance the flexibility of any 
phenomenologically-oriented spectral design. 

3.1 N = 4 example 

A sufficiently representative sample of the spectral patterns may be obtained even via one- 
parametric variations of our matrix boundaries (or, better, row-and-column, vectorial bound¬ 
aries). Some of these patterns may be even revealed at the smallest matrix dimensions IV + 1. 
We found it instructive to use N+l = 5 and to construct Hamiltonians H = H(R) = T+R 
with a single real variable coupling R , with a discrete Laplacean for the kinetic energy (in the 
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notation of Eq. (J2]) we have T = /j( JV+1 )(0)) and with potential 


V {e) 
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where a m = 0 or 1 and T n = 0 or 1. In this way the model is defined in terms of the two 
quintuplets of binary digits forming a complex index 


Q — (<7 4 , u 3 , a 2 , a 2 , <7q) + i {j 4 , t 3 , t 2 , t 3 , tq). 


Let us now pick up one of the simplest indices 


Q a = (0,0, 0,0,1) + i (0,0,0, 0,0,) 

leading to the energy-eigenvalue spectrum {E 0 , E 3 ,, E 4 } defined as roots of secular equation 

E 5 — 4 E 3 + 3E — R = 0. (12) 

This may be re-read as an explicit polynomial definition of the couplings R = R(E) yielding a 
typical wiggly-curve spectral pattern as shown in Fig. [0 The related domain D a of couplings R 
in which the whole spectrum remains real is an interval T> a = (—RH <EP \r^ <EF ^ determined 
by its two KEP boundaries with 

r (kef) = 12 + 8 \/21 y 30 _ 5 ^21 » 1.036340418 . 

125 v 

The ground state remains robustly real so that the latter two KEPs play the role of the points 
of a merger and of a subsequent complexification of the first two excitations E\ and E 2 and of 
the second two excitations E 2 and E 3 at the left and right end of the interval, respectively. 

A different spectral pattern is obtained at 


g b = (1,1,1,1,1) + z (0, 0,0, 0, 0,) 

(see the spectral plot in Fig. [2]). The complexification involves there the ground state E 0 which 
merges with the first excited energy level E\ at the right end of the interval XV Again, just 
the two excited energy levels E 2 and E 3 merge at the left end of physical interval XV 







Figure 1: The R —dependence of energies in potential (TITT) at g = g a . 



Figure 2: The R —dependence of energies in potential (fITT) at g = g 

The domain T> may cease to be compact and/or simply connected. Both of these features 
characterize our third illustration with 


g c = (1,1, 0,0,0) + *(0,0,0, 0,0,) 

for which the Hamiltonian 

^2 + R-l + RO 0 0 

-1 2-10 0 
H C (R ) =0 -12-10 

0 0 -12-1 + R 

0 0 0 -1 2 + R 

yields, for shifted x = E — 2, the factorized secular equation 


(13) 


(x — 1) (x + 1 — R) [x 3 — x 2 R — (3 — R) x + 2 i?] = 0 . (14) 

The spectral pattern becomes anomalous in containing straight lines (cf. Fig. [3]). The coupling 
re-expressed as the function of the energy R(E) develops two poles at x — — 1 and at x — 2. 
This explains why we spot, in the picture, the shared real-constant asymptotics of E 2 (R) and 
E 4 ( R.) at R -C — 1 and of E 0 (R) and E 2 ( R) at R 1, respectively. 
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Figure 3: The R —dependence of energies for Hamiltonian (fT3li at g — g c . 


3.2 N = 6 example 


For the c— subscripted N = 5 Hamiltonian (TT3l) with R — 1 it is straightforward to verify that 
the unavoided energy-level crossing E 2 (R) —> Es(R) —> 3 is accompanied by the parallelization 
of the two related eigenvectors. In other words, the canonical form of the Hamiltonian will 
contain a two-dimensional Jordan block at the R = 1 KEP singularity. From the purely 
phenomenological point of view such a possibility represents an independent and important 
merit of the model which would certainly deserve a deeper systematic study at a general 
dimension N. Before we show in which way such a study gets significantly facilitated by the 
sparse-matrix nature of our present VE— symmetric Hamiltonian (1T0)) . let us still reanalyze 
another, slightly more complicated quantum system described by one of the simplest non¬ 
tridiagonal toy-model Hamiltonians 


H (7 \R) = 


2 -R R - 1 R 0 
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(15) 


This enables us to demonstrate that and how the use of more levels (here, N +1 = 7) and/or of 
more non-vanishing matrix elements in V (here, M( V) = 8) may make the spectrum increas¬ 
ingly complicated and richer (cf. Fig. |4|). 
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Figure 4: The R —dependence of the real eigenvalues of Hamiltonian (fl5lh 

4 Construction of bound states 


4.1 Separable-interaction reinterpretation of V 

In the Dirac’s bra-ket notation let us denote the elements of the orthonormalized basis in 
which we work by ket-symbols |e&) with property V\eo) = |ejv)- This enables us to rewrite our 
VT —symmetric matrix V = H — T of Eq. (HU in the form of a sum 

V = V Z + V ab + V a + Vp (16) 

of four components containing two dyadic terms each, 

V z = -z |e 0 )(e 0 | - z* \e N ){e N \, V ab = a |ejv)(e 0 | + b |e 0 )(ejv|, 

V a = |o)(e 0 | + |ejv)(a| P (N ~ X) , Vp = 0)(e N \ + \e o )0\ V ^ . 

It is easy to verify that our operator V satisfies the PTS alias PHR relation 

y\j)(N+\) _ p(N+ 1) y ' n y) 


What is more important is that the use of projector 

N-l 

n = l f b)( e j 
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enables us to split Schrodinger equation into its two boundary layers and a n—projected sub¬ 
system, 
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We separated the bulk wave function (i.e., its N — 1 inner components |X) = II |^)) from its 
values Xzp = (eo.jvIV’) at the respective spatial boundaries. In terms of an (N — 1)—dimensional 
and, presumably, invertible matrix 

e _1 (f/) = n (ei { n+1) - r (7V+1) ) n (19) 

we may separate the bulk part of our Schrodinger equation, i.e., its N — 1 inner lines 

|a) x- + |/3) x + = T,~ 1 (E) | A") + |ei) x_ + |eiv—1) x+ • 

They may be formally multiplied by the matrix resolvent E(E) from the left yielding the 
II—projection |X) of our wave function in the manifestly E— dependent explicit form 

|X> = E (E) (|a) - | ei )) x_ + E (E) (|/3) - (e^)) x + = \A_(E))x_ + \B + {E))x+ . (20) 

Due to our discrete-Laplacian choice of matrix T the matrix inversion E _1 (i?) —> E (E) is 
non-numerical (see the next subsection). Hence, vector (I2U1) is known in closed form and it 
may be inserted in the remaining two lines of Schrodinger equation. We end up with the two 
homogeneous linear equations 

{2-E-z)x. + {l3\V {N - 1) \X)-{e l \X) + bx + =0, (21) 

a x. + (a\V [N ~ l) \X) - (ejv-i \X) + (2 - E - z*) x+ = 0 (22) 

determining the last two missing components x± of the wave function. Naturally, the set of 
bound state energies E is finally obtained as zeros of the related secular two-by-two determinant. 


4.2 Non-numerical diagonalization and construction of E(£ l ) 


Undoubtedly, the exact solvability of our present model is due to the non-numerical invertibility 
of the (M + 1) by (M + 1) matrices (TT9l) at any dimension M + 1. More details may be found 
elsewhere [21]. For our present purposes it is sufficient to remind the readers that the Chebyshev 
polynomials U n (y) of the second kind satisfy the relation [22] 
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provided only that we choose the value of y as one of the M + l roots of the polynomial equation 

Um+i(v) = Le -> 


(M+i) Q (.3 + 1)tt 

y = y) = 2 cos ■ 




(24) 


M + 2 ’ 

As long as all of these roots are known in closed form, we may arrange the (M + 1)—plet of 
the column-vector solutions of Eq. (I23l) into a unitary matrix U and diagonalize the discrete- 
Laplacean matrix 7d M+ b since 
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(25) 


with 5j R1+1) = 2 ^1 — Uj M+1 ^ • Once we now choose M = N — 2, it becomes easy to return to 
Eq. (jT9|) and to perform, non-numerically, the inversion E -1 (E) —y E (E). 


5 Construction of physical metrics 

Several detailed studies [23] confirmed, by non-numerical constructive means, that in an ad hoc 
physical Hilbert space the evolution of the differential-equation model (EJ) + (EJ) can be 
made unitary. In Ref. [8], analogous results were obtained for the difference-equation model 
dSJ) + © ■ hr both of the latter contexts the explicit construction of the physical Hilbert space 
i.e., of its inner-product metric operator 0 represented not only the quantum-theoretical 
necessity but also by far the most difficult mathematical challenge. 

5.1 Expansions of © in terms of eigenvectors of iC 

In technical terms one has to find such an operator 0 = 0^ > 0 which is connected with the 
given Hamiltonian H ^ H t via the S'—space Hermiticity condition 

#10 = 0 # ( 26 ) 

In an abstract context of pure mathematics such an S'—space Hermiticity relation was studied 
by Dieudonne m- More or less simultaneously the same generalization of Hermiticity proved 
useful for physicists who found a number of its useful practical applications in many-body 
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quantum systems mm- In what follows we shall show that the construction of a suitable 
Hermitizing metric 0 remains feasible for our present Hamiltonians (fTOlh 

For finite-dimensional Hamiltonians there exist two alternative and efficient methods of the 
reconstruction of the necessary Hermitizing metric 0 = O(ih). The more universal spectral- 
expansion approach as discussed in Ref. [26] is based on formula 

N 

Q = Yl \ n )) K l(( n \ , KnGM\{0}. (27) 

0 

The special kets |n)) are defined here as a complete set of solutions of the conjugate-Hamiltonian 
Schrodinger equation 

H ] | n)) = E n \n )), n = 0,1,..., N . (28) 

In general, one merely has to employ the method of solution as described in section [I] above. 
Alternative versions of the method may be developed for some simpler special cases. One of 
them will be described in what follows. 


—* —> 

5.2 Closed-form metrics for \/3) = 0 

Let us restrict our attention to the subset of models m with trivial \/3) = 0, i.e., with the 
conjugate Schrodinger operator 


(h^ n+1) - E I^ N+1 ^ 
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where we abbreviated 2 — E = 2y and partitioned 




(u(y) | 

a 

. -s-'iv) 

l_ 


(29) 


(30) 


We insert this matrix in Eq. (T28|) and get a re-partitioned conjugate Schrodiner equation 


{u(y n ) I a 

~S(y n ) V^\u(y n )) 


n)) = 0, 



(31) 
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As long as matrices H and W are isospectral, the bound-state energy E n (together with its two 
reparametrizations 2 y n = 2 — 2 E n — t n + l/t n ) is already known, at any n — 0,1,. .., N, from 
section Q] Thus, the hrst line of Eq. (T3T1) may be omitted as redundant. In terms of an arbitrary 
normalization constant x n the rest of this equation determines the ultimate N— dimensional 
vector 

\Xn) = S~ l {y n ) V {N) \ u(y n )} x n . (32) 

It is important that the explicit form of the inverse of our triangular N by N auxiliary matrix 
is obtainable in closed form, 


S~ 1 {y n ) = U{t n )U(l/t n ), U{r) 


1 r t 2 

0 1 r 

0 0 1 


.N -1 " 
.N— 


: ■ ■ ■ ■ ■ ■ t 

0 0 ... 0 1 


(33) 


Thus, the insertion of Eqs. (l32l) and (133|) in Eq. (1271) completes the construction of the metric. 
We see that due to the simplification | (3) =0 we did not have to pre-diagonalize and invert the 
tridiagonal kinetic-energy matrix T at all. 


6 Sparse-matrix metrics 

The key advantage of the strategy of the preceding section may be seen in an immanent 
guarantee of the positivity of the metric. A weak point of the method is that even if the 
Hamiltonian itself admits an exceptionally simple (i.e., diagonal or other sparse-matrix) form 
of the metric 0 = Q(H) ^ /, it is difficult to identify such a most welcome solution of 
Dieudonne’s Eq. (1261) as a special case of the general spectral expansion (T27|) . A different 
construction strategy is needed in such a situation. 

6.1 Recurrent constructions 

In the spectral-expansion formula HTTP the physics-determining choice (i.e., the S— superscripted- 
Hilbert-space-determining choice) of the (N + 1)—plet of positive parameters is fully at the 
user’s disposal (cf. Ref. [IT] for an exhaustive analysis and explanation and physical interpreta¬ 
tion of such an apparent ambiguity). Naturally, in an alternative approach one may choose free 
parameters directly as some matrix elements of 0. Then, one has to determine the remaining 
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elements of the Hermitizing metric 0 via a direct, brute-force solution of the linear algebraic 
system of (N + l) 2 equations (1261) (which are not always mutually independent). 

An efficient implementation of the latter recipe is much less universal and its feasibility 
depends strongly on the structure and properties of the input Hamiltonian H ^ HT For 
illustration purposes let us consider, therefore, just the one-parametric five-dimensional real- 
matrix Hamiltonian H C (R) of Eq. (fl3jl . In accord with Eq. (TT4| and Fig. [3] we know that the 
energies are all real, say, at any R < 1. It is also easy to prove that the simplest possible 
Hermitizing metric @( c ) > 0 may be chosen as a diagonal matrix with elements 0^ = 1, 
= 022 = 033 = 1 — R and ©44 = (1 — R) 2 . Nevertheless, it would be fairly difficult to 
derive this result from the mere inspection of spectral expansion f)27)l . An alternative, recurrent 
recipe may be recommended. 

Lemma 1 The complete family of the real and symmetric matrices 0 which Hermitize Hamil¬ 
tonian H C (R) of Eq. (Q3j) may be constructed by solving Eq. (d6j) in recurrent manner. 

Proof. 

Let us treat matrix elements u = 0 O i, z = © 02 , q = @ 03 , P — ©04 and t = @22 as a quintuplet 
of independently variable lower-case parameters in the metric. The subsequent inspection of 
relations fl26|h viz., 

4 4 

■M-ij = ^ ik _ ®inH n j = 0 

k =0 n =0 

reveals that their ^ = 04 item defines ©i 4 = q — Rq, their j j = 03 item defines ©i 3 = R q + z + p 
and their ^ = 02 item defines @i 2 = Rz + u + q. Next, after insertions of these elements their 
ij = 14 item defines 0 2 4 = z — Rz and their ^ = i3 item defines 0 23 = q + Rz + u. Next, after 
insertions of these elements their = 24 item dehnes ©34 = u — Ru. Next, after the insertion 
of this element their ^ = 12 item defines diagonal element ©n = — Rq — p + t — Rz and their 
ij — 23 item dehnes diagonal element © 33 = — Rq — p + t — Rz. Finally, after insertions of 
these elements their VJ = 0 i item dehnes ©oo = (R z — t + z + Ru + R q + p)/(—1 + R) and their 
ij = 34 item dehnes ©44 = —z + Rp — Rq + R 2 q — p — Ru + t — t R + R 2 z + R 2 u. 

□ 

In Fig. |3] we may notice (and, if asked for, we may readily prove) that for Hamiltonians 
H C (R ) the energy spectrum is real in the interval of R. G (— 00 , R ma x ) with boundary value 
Rmax ~ 1.065260704 which is larger than one. At the same time, a return to the proof of the 
preceding Lemma implies that every candidate 0 = 0{t,u,z,q,p ) for the metric (say, with 
t — 1) is a diagonal matrix if and only if one sets u = z = q = p = 0. The resulting matrix 
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0(1, 0,0, 0,0) is positive definite (i.e., it may serve as a physical metric) if and only if R < 1. 
This is an incompatibility which has the following important physical interpretation. 

Lemma 2 If we assign the same diagonal metric 0(1, 0,0, 0,0) to Hamiltonian H C (R) at all 
R < l, the passage through R = 1 will necessarily lead to a phase transition, not accompanied 
by the loss of observability of the bound-state energy spectrum. 

Proof. 

Firstly, it is obvious that the value of R = 1 is a KEP singularity of matrix H C (R). At this 
point, secondly, the diagonal physical metric ceases to be positive definite. We may conclude 
that for R G (1 ,R max ), any acceptable metric 0 = Q(t,u, z,q,p) becomes necessarily non¬ 
diagonal. To the left and to the right from the KEP singularity R — 1, the metrics 0 (i.e., 
the physical Hilbert spaces PL^) are necessarily different so that also the quantum system in 
question must be necessarily characterized by the different complete sets of the operators of 
observables. 

□ 

At the right end of the smaller physical interval of R G (1, R m ax ) the mathematical as well 
as physical nature of the KEP singularity is less unusual since at this value the two highest 
energies cease to be real. They merge and, subsequently, complexify. The merger value of 
energy E = x + 2 may be determined as a root of quartic polynomial x 4 — 2 x 3 — 3 x 2 + 6 = 0. 
Numerically one reveals that E 3 —> 3.233152848 <— E± in the limit of R —> R max - 

6.2 Diagonal metric candidates 

The recurrent construction method as sampled in preceding subsection leads to the metric- 
operator candidates 0 in which certain merits (e.g., sparsity) are combined with shortcomings 
(e.g., a possible premature loss of the necessary positivity). Naturally, at any dimension N + 1 
the merits get maximal and, simultaneously, the shortcomings are minimal for diagonal matrices 

0 = Qidiag) 

so their subset deserves an enhanced attention. 

At any matrix dimension, the same partitioning may be applied to Hamiltonians H^ N+ 1 ) = 
T( n+ i) + V^ N+l \z,a,b,a, fd) of Eq. (fTOjl as well as to the related diagonal metric candidates. 
Once we separate the outer elements we may abbreviate @qo“^ = m 0 and and 

use the subscripts j = 1,2,..., N —1 numbering just the inner part of the diagonal = de¬ 

forming a smaller, tilded diagonal matrix d. 

We shall assume that the Hamiltonians as well as the metric candidates are real matrices. 
In the light of the IV = 4 results of the preceding subsection we feel encouraged to classify, at 
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all TV, the scenario in which the two matrices remain compatible with the hidden Hermiticity 
constraint (l26j) . 

First of all let us restrict our attention to the inner, (TV — 1) by (TV — 1) matrix partition of 
Eq. m- The potential cannot contribute so we must set d\ — d 2 , ■ ■ ., d^~ 2 = d^-i- Without 
any loss of generality we have dj = d > 0 at all j, i.e., d = dl 

After we turn attention to the four outer partitions of Eq. (I26jh we may study just its first 
and last columns of the real and symmetric matrix set A4 = //'(-) — 0 II = 0 with a trivial 
main diagonal. The Erst column offers a set of TV relations 

(/3at_i — 1) m 0 = d (a 1 — 1), /3 N ~j m 0 — daj, j = 2,3,..., TV — 1, b m 0 = m N a . (34) 

They have to be complemented by the TV—plet of the last-column relations 

arriN = m 0 b , Q>N-jrriN = df3j , j = 1,2,..., TV — 2 , (an — 1) tum = d (/3?v-i — 1) (35) 

in which, due to the symmetry of A4, the Erst item is already redundant of course. 

As long as we must have positive d > 0 and > 0 for any diagonal metric, the former 
set of relations may be now read as the definition of the TV—plet of “acceptable” parameters 
atj and a in terms of the respective “optional” parameters @N-j and b. I 11 such a case we may 
insert these quantities in the latter set of relations and discover that up to the uninteresting 
trivial case (with a completely vanishing nonlocal-boundary-condition part of the Hamiltonian) 
the whole resulting set of relations degenerates to the single sufficient condition moupv = d 2 . 
This completes the proof of the following result. 

Proposition 3 For the (TV + l) by (TV + 1) real-matrix Hamiltonian HTU\) there exists a positive- 
definite diagonal metric if and only if the N—plets of the dynamical parameters a (and 

a) and (3 (and b) satisfy the N—plet of linear relations |^| ) where d > 0 is a fixed inessential 
overall multiplier and where momjv = d 2 . 

Remark 1 While the diagonal metric may vary with just a single “essential” free parameter, 
the variability of the energies remains controlled by an N—plet of free couplings. 

7 Parametric domains J of positivity of the metric 

7.1 TV = 4 example 

A decisive merit of the above-sampled and comparatively complicated recurrent construction 
of the complete set of the metric candidates should be seen in the facilitated construction of 
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their sparse-matrix subsets and even, in some cases, of their (2 K + 1)—diagonal band-matrix 
versions. This is well illustrated by the above-derived explicit formula 


m(00) 

u 

z 


p 

u 

m(ll) 

Rz + u + q 

Rq + z + p 

q- Rq 

z 

Rz + u + q 

t 

Rz + u + q 

z — Rz 

Q 

Rq + z + p 

Rz + u + q 

m( 33) 

u — Ru 

p 

q — Rq 

z — Rz 

u — Ru 

m( 44) 


for H = H C (R) of Eq. (fT3l) where we abbreviated 
Rz-t + z + Ru + Rq + p 

m(00) = -—-, m(ll) = m( 33) = —Rq — p + t — Rz 

—1 + R 


and 

m(44) = — z + Rp — Rq + R 2 q — p — Ru + t — t R + R 2 z + R 2 u . 


(36) 


(37) 


Almost all of the matrix elements are just linear functions of the parameters t, u, z, q,p and R. 
For this reason, the sufficiently small deviations of the parameters from their diagonal-metric 
limit u = z = q = p = 0 can merely restrict (but not entirely destroy) the robust positivity 
(i.e., acceptability) of the diagonal-metric candidate. 

For illustrative purposes let us pay attention just to the first nontrivial, tridiagonal-matrix 
special case of Eq. ()36jl and let us re-analyze the one-parametric family of the candidates 
0 = 0(1 — R, u, 0, 0, 0) for the metric. Our Erst observation is that even at the very small 
fixed u we cannot exclude irregularities in the R —> 1 limit. The reason is that at any fixed 
value of u > 0 the matrix element m(00) as given by Eq. (15T|) remains singular in the limit, 
lim^!- m(00) = — oo. Near the singular value of R the matrix element m(00) is dominant 
and becomes approximatively equal to one of the eigenvalues. Simultaneously, as long as this 
matrix element becomes large and negative for positive u in this limit, the candidate matrix 
becomes indeterminate and, hence, inaeceptable near the KEP boundary of R ~ 1. 

Fortunately, after one regularizes the R — 1 singularity in m(00), the rest of the metric 
candidate already appears to be a smooth function of its parameters. At their small values 
one may expect just small deviations of eigenvalues from their unperturbed limit with 9 0 = 1, 
with the degenerate triplet of the linearly R— dependent eigenvalues 9 1 >2 ,3 = 1 — R and with 
the quadratically A—dependent eigenvalue 6*4 = (1 — R) 2 . 

Once we consider just the Erst nontrivial and properly regularized tridiagonal matrix 0 = 
0(1 — A, (1 — R) w, 0, 0, 0), our expectation of the smoothness of the variation of the eigenvalues 
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Figure 5: The R —dependence of the eigenvalues of the tridiagonal N — 5 matrix @(f, u, 0, 0, 0) 
of Eq. (1361) at normalization t — 1 — R and small u — (1 — R)/ 10. 



Figure 6: The same R —dependence as in Fig. [5] but at almost critical u = 6(1 — R)/ 10. 

are conhrmed. This is documented by Fig. fusing a small w = 1/10. We see there that 
whenever R < 1 the spectrum of the metric remains positive. As expected, its R —dependence 
keeps the clear trace of the separation of one eigenvalue of an almost constant type from an 
almost linear triplet and from an almost quadratic shape, at the not too small values of R at 
least. Only far to the left from R — 1 (i.e., near R — 0) this pattern changes since avoided 
crossings enter the scene. 



Figure 7: The u —dependence of the eigenvalues of 0(1 — R , u, 0, 0, 0) at R = 8/10. 

With the growth of w towards a certain not too small w = w cr iti ca i we expect that at least 
one of the eigenvalues of the tridiagonal metric candidate 0(1 — R,(l — R) w, 0, 0, 0) will vanish 
so that beyond the boundary point w cr i t i ca i the candidate matrix will become indeterminate 
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(i.e., for our present purposes, useless). Recalling and comparing Figs. [5] and [6] we may visualize 
the empirical observation that the point R m .i n of the loss of the positivity of the metric candidate 
moves up with the growth of w. 

Using the brute-force numerics we managed to arrive at a rough estimate of the maximal 
possible w = w cr iticai £ (7/10,8/10). Nevertheless, after a return to Fig. [6] one sees that in a 
comparatively large vicinity of w cr iti ca i the minimal positive eigenvalue 9 o is in fact rather small 
in the whole admissibility interval of R £ (R mm , 1). Thus, the numerical invertibility of the 
metric would certainly deteriorate long before the value of w reaches its critical value. Thus, 
the model will only prove numerically well conditioned well below the maximal w cr iu C ai- 

The regularization u —> w is not needed at a fixed dynamical coupling R < 1. We choose 
R = 4/5 and displayed, in Fig. [TJ the graphical proof of the positivity of all of the eigenvalues 
{9j} of the tridiagonal metric candidate 0(1/5, u, 0, 0, 0) in an interval of the length which is 
not too small (i.e., in fact, it is bigger than 1/5). 

7.2 Arbitrary dimensions N + 1 

In Ref. [27] it has been shown that the recurrent constructions of the metric candidates remain 
particularly efficient for the sparse-matrix Hamiltonians which are at most tridiagonal. Hence, 
in a search for a more challenging/rewarding illustrative example let us now return, once more, 
to the asymmetric and non-tridiagonal real-matrix Hamiltonian H^ 7 fR) of Eq. (1I5|) . For this 
model we know that the spectrum is real in Hermitian limit R —)■ 0 [21] and that it remains all 
real at small R ;$ 0.555 at least (cf. Fig. [4] above). This encourages us to formulate and prove 
the following generalization of lemma [H 

Proposition 4 The complete family of the real and symmetric matrices @( Ar+1 ) which pseudo - 
Hermitize the minimally nonlocal-boundary M(V ) = 8 Hamiltonians H^ N+1 1 as sampled by 
matrix H^ 7 \R) of Eq. / f731) may be constructed by solving pseudo-Hermiticity condition [26 1) in 
recurrent manner. 

Proof. 

Proof proceeds in a complete parallel to the proof of lemma [TJ It starts by the interpretation 
of the Erst row of matrix elements aq = 0 O o, x 2 = ©cm • • ■, aqv+i = ©o n ( one has, in our new 
illustrative example, N = 6) as independent parameters. After the explicit display of relations 
(1261) in matrix form 

N N 

M t , = y e kj - y e tn H nj = 0 

fc =0 n =0 
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one has to proceed as above. Thus, the selection of the set of independent equations starts 
from M.qn = 0 and moves to the left up to the diagonal element which is found trivial and 
omitted. In the next round one moves down and repeats the selection. In this manner one 
obtains the sequence of (N + l)N/2 linear equations Aiij = 0 (with i < j) for single unknown 
quantity yielding easily, with insertions performed after completion of each row, the sequence 
of the required explicit formulae for the respective matrix elements ©*+ ij(xi,X 2 , ■ ■ ■, Xjv+i)- 

□ 

As a result of the application of the above algorithm one has the complete, multiparametric 
set of matrix candidates 0( JV + 1 )(o;i, x 2 , • • •, apv+i) for the metric. All is now prepared for the 
final step of the construction of the metric. It is sufficient just to specify a non-empty domain 
J of parameters Xi,x 2 , ■ ■ . ,£jv+i (be., of vectors x) in which the acceptability (i.e., positivity) 
of the candidate matrix 0 is guaranteed. 

At all N and at all of our present nonlocal boundary conditions the latter task is facilitated 
by the fact that at z = 0 and in the Hermitian limit of a —> 0 and /3 —>■ 0 the Hamiltonians 
degenerate to the discrete Laplacean T (Ar+ b with real and non-degenerate spectrum. This 
means that we may choose a trivial metric 0 = / in this limit. Next, in a sufficiently small 
vicinity of the values of 0 = 0, of a = 0 and of f3 = 0 we are sure that the energy spectrum 
of H remains, for continuity reasons, real and discrete (indeed, the complexihcation may only 
proceed via a smoothly proceeding merger of at least one energy doublet [TO]). By construction 
and for similar reasons, also the Hermitizing metrics (with elements being linear functions of 
parameters x) will deviate from the unit matrix while still remaining positive definite in a 
suitable non-empty domain J. We just completed the proof of the following result. 

Theorem 1 At the sufficiently small dynamical parameters z, a and (3 the metric candidates 
0 (A r +i)(£) 0 f p r0 p 0S iti 0n ^ constructed as perturbations o/ 0 (JV+ b(0 ) = / will be positive def¬ 
inite (and will Hermitize the Hamiltonians) in a non-empty domain J of sufficiently small 
parameters x. 

In order to illustrate the latter result let us pick up again Hamiltonian H^fR) of Eq. (fT5]i 
and using the algorithm of Proposition [T] let us construct the following illustrative sparse-matrix 


22 



Figure 8: The R —dependence of the eigenvalues of metric candidate 
candidate for the metric, 


e (7) (l,0,0,0,0,0,0) = 
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(38) 


The variations of its eigenvalues with R are displayed in Fig. [8] One may conclude that besides 
the required Hermitian-limit property liniR^o @^(1, 0, 0, 0, 0, 0, 0) = I* 7 - 1 , the matrix remains 
also safely positive definite in a not too small vicinity of the Hermitian-Hamiltonian limiting 
point R = 0. 

Another, different lesson taught by the latter example is that at any dimension N, the 
breakdown of the tridiagonality of the Hamiltonian may make the structure of the metrics less 
predictable. In a way guided by our N = 6 illustrative O jk one may expect, for example, 
that even in the maximally sparse cases the metrics may not possess the few-diagonal matrix 
structure anymore. 


8 Pseudometrics 

The pseudo-Hermiticity constraint 

H ] V = VH (39) 

is formally equivalent to the hidden Hermiticity constraint (126)1 . with two differences. Firstly, 
the matrix of metric 0 in Eq. (126|) is always (perhaps, just tacitly) assumed positive definite 
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and, secondly, it is allowed to be Hamiltonian-dependent. In contrast, the (invertible) operator 
of pseudometric V in Eq. (j39|) is usually required Hamiltonian-independent and indefinite (plus, 
sometimes, such that V 2 = I). 

Without taking these differences into account, it is obvious that among the sparse-matrix 
metric candidates one should find all of the indefinite pseudometrics. In particular, using the 
results of preceding paragraph 17.21 on the non-tridiagonal sample Hamiltonian (TT5|l it is easy to 
verify that during an exhaustive algebraic solution of Eq. (126|) one obtains the result 

0 (7) (0,0, 0,0, 0,0,1) = V (7 \ 

i.e., one really re-obtains also the initial Hamiltonian-independent “input” matrix ([2]) of parity. 

In a more general methodical context it seems tempting to weaken the constraints im¬ 
posed upon the pseudometrics. In the PHR framework of review [10] precisely this type of 
generalization (i.e., a replacement V —> V such that not necessarily V 2 = /) was proposed 
by Mostafazadeh. For our concrete toy-model H^\R) of Eq. (fT5]h for example, this would 
mean that whenever necessary, one could search for several alternative, multiparametric and 
Hamiltonian-dependent sparse-matrix pseudometrics. 



Figure 9: The R —dependence of the eigenvalues of the pseudometric candidate fjTOjl . 

Obviously, the most natural pseudometric candidate would be matrix ©T)(0, 0, 0, 0, 0,1, 0) 
which has the compact form 
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Figure [9] demonstrates that at non-vanishing and not too large R / 0 this matrix is invertible 
and, hence, it could play the role of an indefinite pseudometric. Still, a serious practical 
drawback of the latter pseudometric candidate with eigenvalues r lies in its ill-conditioned 
invertibility caused by the smallness of the middle eigenvalue 73 in a fairly broad vicinity of 
R — 0. For a remedy we may replace Eq. (140]) . say, by the following one-parametric family of 
more general candidates for the Krein-space pseudometric, 

V = V{Ci = © (7) (0,0,0,0,0,1,0) ± £R . 

We could search for a value of £optimal which would (approximatively) maximize the distance 
of the smallest eigenvalue of such a generalized pseudometric candidate "P(£) from zero. In a 
preparatory step we shall recall the approximate R —independence of r 3 in Fig. [9]at small R and 
at £ = 0. This enables us to simplify our optimization problem by replacing variable R by zero. 
In other words, a good guess of £optimal will be provided by the value of £approximate obtained 
after the replacement of the true secular polynomial /(£, R, £) = det (v(£) — ( J (jV+1 ^ by its 
reduced two-variable form /(£, 0 , r). 

A serendipitious benefit of the R — 0 simplification arises via the following closed-form 
factorization formula 

/«,0,T) = (T + £)((T + f) 2 —2) (((r —f) 2 —2) 2 — 2 ) . (41) 

This formula enables us to find the recommended value of £ (which would give all eigenvalue 
roots r sufficiently far from zero) in two alternative and in absolute value equally large closed 
forms £ approximate = ±1/a/4 + \/8 ~ ±0.38 of the point of intersection of the fourth root 
t 3 (£) = —£ of /(£, 0,t) with the third root 72 (£) = £ — a/2 — y/2 or with the fifth root 
t 4 (£) = £ ± a/ 2 — \/2, respectively. 


9 Summary 

In the conclusion one may predict that our present success in keeping the multi-parametric 
discrete square well still friendly and tractable by non-numerical (or, at worst, semi-numerical) 
means could have opened a way towards its further future study. Not only in the obvious 
mathematical direction towards the continuous-limit bound-state regime with N = oo but also 
towards establishing its phenomenological relevance. In the dynamical context of scattering, 
after all, the first results and an independent inspiration could be already found in Ref. [15] . 

Another, independent, implementation-oriented appeal of our present general endpoint in¬ 
teractions may be found in the context of the well known phenomenological Su-Schrieffer-Heeger 
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Hamiltonians, the tridiagonality of which is often strongly violated by an addition of elements 
in the upper-right corner and/or in the lower-left corner of the Hamiltonian matrix (cf., e.g., 
Ref. [28]). Naturally, such an exceptionally remote endpoint interaction already acquires a new 
physical meaning. Indeed, the main purpose of its absence or presence gets shifted from its 
present act ion-at-a-distance interpretation to the possibility of control of the transition between 
the open-chain model of Eq. ([T]) and its topologically nonequivalent closed-loop alternative. 

The simplest sample of the present endpoint interactions in such a new, spatial-topology- 
related role would be provided by the use of the mere two extreme elements b N ^ 0 and d 0 ^ 0 
in interaction matrix Eq. (J8]), or of a ^ 0 and b ^ 0 in Eq. (TTOl) . The addition of more (or, 
perhaps, all) non-vanishing endpoint elements would lead, in this sense, to a simulation of a 
toroidal geometry of the interactions. An even stronger emphasis put on the similar simulation 
of topological effects via long-ranged alias endpoint-matching interactions may be found in the 
context of quantum graphs [29]. One can summarize that in all of the similar multiparametric 
families of mutually related models the information about the change of the topology and/or 
spectrum is carried strictly by the nonlocal boundary conditions. 

Needless to add that in many physical applications of the THS pattern (like, e.g., in 
the so called interacting boson models of heavy atomic nuclei DU or in the descriptions of 
the schematic non-Hermitian Lipkin-Meshkov-Glick manybody model [TT[ [30], etc) the lo- 
cality/nonlocality is not an issue at all. After development and description of methodically 
oriented nonlocal toy models, one may feel encouraged to replace the F —local external forces 
V (£) by their various F— nonlocal, integral-kernel generalizations R(£, £') also in some more 
realistic toy-model considerations. In the related future work it will only be necessary that the 
operators of observables still remain sufficiently user-friendly. 
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Appendix A. Operators of observables in non-Hermitian 
representations 

According to the conventional model-building strategy the physical Hilbert space W p> of quan¬ 
tum theory is chosen, a priori , in one of its simplest realizations. Typically, one decides to work 
with the space L 2 (M) of all square-integrable complex functions of coordinate or momentum. 
On this background, the Bender’s and Boettcher’s discovery [31] of the possibility of a peaceful 
coexistence of observability with non-Hermiticity changed our perception of what is an optimal 
mathematical description of a quantum system. 

The deeper and more detailed account of the amended theory and/or of its various innovative 
applications may be found elsewhere mm- It has been clarified that a (mathematically more 
or less trivial) replacement of the “primary” Hilbert space 'H (p> by its “secondary”, unitarily 
equivalent alternative "HA) might exhibit multiple specific merits even if the latter physical 
Hilbert-space requires a nontrivial definition of the inner product and even if it ceases to 
belong, therefore, among the simplest possible realizations of the picture. 

In the new theory one has to proceed, basically, in two steps. Firstly, one has to replace 
the initial (but, by assumption, prohibitively complicated) physical Hilbert space by 

its intermediate substitute l-lS F > which is manifestly unphysical but, by assumption, much 
friendlier. In the notation as introduced in [12] this is the first step towards simplification 
which is to be achieved by means of a suitable invertible mapping fh 

In the second step one takes into account that whenever the auxiliary map itself appears 
non-unitary, the new representation H of a hypothetical Hamiltonian f) which was selfadjoint 
in the prohibitively complicated becomes non-Hermitian in unphysical 

H = ^ fiff) (fi -1 )* = H ] . (42) 

From the latter relations we may eliminate f), introduce an abbreviation = 0 and get the 
original Hermiticity relation f] = rewritten in the form of Dieudonne’s constraint (l26|h Thus, 
in many cases of interest we may define a new operator H * = Q~ 1 H^ 0, require that H = 
and declare our new Hamiltonian self-adjoint in the third Hilbert space endowed with the 
ad hoc inner-product metric operator 0 / /. 

The S'—superscripted space coincides with ji^ F > as a topological vector space. These two 
Hilbert spaces merely differ from each other by the different definitions of the respective Her- 
mitian conjugations, i.e., of the respective inner products, 

IW'iW'2)] <F, #[W'iIV- 2)] <S) = [«’i|0|>fe}r ) • (43) 
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In opposite direction, once we fix the metric (with properties 0 = @f > 0 etc [fl]), any 
candidate A for an operator of an observable must be essentially self-adjoint in This 

means that we must demand either that A = A* (when considered in the S —superscripted 
physical Hilbert space), or that 

A t 0 = 0A (44) 

when our quantum system is represented in the F— superscripted, calculation-friendly space 
(which is the preferred choice, after all, in the vastest majority of applications). 


Appendix B. An explanatory remark on the concept of 
locality in non-Hermitian context 

The current acceptance of the compatibility of non-Hermiticity with observability was perceiv- 
ably accelerated by the Bender’s and Boetcher’s pH] restriction of attention to Hamiltonians 
of ordinary differential PTS form 

H {BB) = -^ + v(BB) (0- (45) 

d £ 2 

These Hamiltonians were defined in the most common Hilbert space L 2 (M), but this space (i.e., 
in our present notation, had to be reclassified as unphysical. One must Hermitize H ( ' BB ' ) 

in by using an amended metric © 7 -I. 

One of the most important phenomenological consequences of the acceptance of the PTS 
theory is that in operators (1451) as well as in wave functions ^>(£), the variable £ does not carry, 
in general, the habitual meaning of an observable point-particle coordinate anymore (cf. [32]). 
The reason lies in the generic violation of its observability status as given by condition f|44|) . 
Indeed, for any operator of position Q defined in we may immediately rewrite the latter 

condition as follows, 

Q = fT 1 qfi^Q t = fi t q(fr 1 ) t . (46) 

In the next step we turn attention to space in which we are allowed to use the principle 
of correspondence. This enables us to identify the self-adjoint operator q with an operator of 
position. I 11 the final step we notice that its spectrum (i.e., real eigenvalues q) coincides with 
the spectrum of Q so that up to the trivial cases characterized by the diagonal-matrix Dyson’s 
maps the property of locality of an operator of interaction in 'H (P) will translate into its 
nonlocality in and vice versa. 



In a constructive manner the latter apparent paradox has been clarified in Ref. [20] where 
the authors emphasized that the concept of a point-particle position alias spatial coordinate can 
only be consistently introduced in the physical Hilbert space with trivial metric (i.e., in 
In opposite direction the illustration of the paradox may be provided by the most popular 
Bender’s and Boettcher’s toy-model (1451) in which the potential V was chosen local in 
unphysical space Still, the model itself remains nonlocal in both of the physical spaces 

(where the nonlocality is due to the non-diagonality of in Eq. (|46|P and (where 
the nonlocality is due to the non-diagonality of product 0 = fhQ [32])■ We may summarize 
that in the model-building context, the locality/nonlocality of non-Hermitian V in does 
not change the generic nonlocal nature of the physical lower-case interaction operator 


nonlocal ) _ q y (local/nonlocal) 


-1 


(47) 


which acts in 'H [P) and is Hermitian there. Once deduced from the PTS models, this interaction 
term is generically nonlocal, i.e., almost always nonlocal, up to a few very rare exceptions like 
the one described in Ref. [33]. 

Naturally, the generic nonlocality of the action of forces (147|) is a serious obstacle for their 
being put in existence in experimental setups. Even for theorists, such a spatial nonlocality of 
action appears counterintuitive j!8j . mainly due to the widespread habit of keeping the contact 
with point particles (and, hence, of using the principle of correspondence) - which need not 
remain a successful strategy even in quantum models which remain safely Hermitian [31]. 

Incidentally, the clear-cut separation of the physical concept of observable locality (taking 
place in from its purely formal parallel in offers also an answer to the authors of 

the recent letter [35] who claimed, erroneously, that their results “essentially kill any hope of 
PTS quantum theory as a fundamental theory of nature”. The core of the misunderstanding 
clearly lies in the fact that the authors of loc. cit. deduced their message from an extremely 
unfortunate mix-up of the concepts of locality in T-i (P \ in and in u {F) [3S3- 

Let us add that the same inessential role was also played by the locality of non-Hermitian 
potentials V(£) in the effective, incomplete-information quantum models of Refs. [18], T9] in 
which a non-unitarity of the scattering has been admitted. In our present, more conventional, 
unitary bound state setting the key benefit of the above analysis lies in the conclusion that the 
locality of potentials V^ BB \tf) did not have any immediate observational meaning and it was 
only relevant due to the facilitated solvability of the related Schrodinger equation. 
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